Ultrafast Optical Excitation of a Persistent Surface-State Population in the 

Topological Insulator Bi2Se3 
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Using femtosecond time- and angle- resolved photoemission spectroscopy, we investigated the 
nonequilibrium dynamics of the topological insulator Bi2Se3. We studied p-type Bi2Se3, in which 
the metallic Dirac surface state and bulk conduction bands are unoccupied. Optical excitation 
leads to a met a- stable population at the bulk conduction band edge, which feeds a nonequilibrium 
population of the surface state persisting for >10 ps. This unusually long-lived population of a 
metallic Dirac surface state with spin texture may present a channel in which to drive transient 
spin-polarized currents. 

PACS numbers: 78.47.J-, 73.20.-r, 72.25.Fe, 79.60.-i 



Three dimensional topological insulators (TI) are fasci- 
nating compounds that combine an insulating bulk elec- 
tronic band structure with a surface state (SS) that 
crosses the Fermi level Ep [1]. The SS has novel proper- 
ties such as topological protection \% |3| , suppression of 
backscattering [4] , and a helical spin texture [5|, l6| . Its re- 
sponse to optical excitation is expected to lead to several 
interesting phenomena, such as large Kerr rotation [3|, [D] 
and surface spin current [8[. However, few experimental 
studies have focused on the SS response due to the fact 
that surface signatures are easily overwhelmed by bulk 
contributions to transport [9|, [lOl] and optical measure- 
ments [ll|,[l2|l. 

In this Letter we report time- and angle- resolved pho- 
toemission spectroscopy (trARPES) experiments on the 
topological insulator Bi2Se3 and elucidate the electron 
dynamics in response to optical excitation 14l4l6j. The 
femtosecond time, energy, and momentum resolution al- 
lows us to probe fundamental scattering processes di- 
rectly within the electronic band structure, and thus dis- 
tinguish the separate dynamics of the bulk and surface. 
We use p-doped samples to start with a completely unoc- 
cupied SS. Optical excitation populates high- lying bulk 
states, which rapidly decay to lower energy states via 
inter- and intra-band phonon-mediated scattering pro- 
cesses. Within 2 ps a meta-stable population forms at the 
bulk conduction band (BCB) edge. This population acts 
as an electron reservoir which fills the SS with a steady 
supply of carriers for ~10 ps. The persistent occupation 
of a metallic state is a unique situation, since the absence 
of a bandgap in metals typically means that there is no 
barrier to rapid recombination. Further, the SS has the 
same Dirac dispersion seen in equilibrium experiments, 
where a novel spin-momentum locking has been demon- 
strated [6] . We propose that this long-lived spin-textured 



metallic population presents a channel in which to drive 
transient spin-polarized currents. 

Our trARPES setup is described in detail in the sup- 
plemental material [13]. Single crystals of Bi2Se3 were 
synthesized by slow cooling a binary melt, with nomi- 
nal Bi:Se ratio of 35:65. p-type Bi2Se3 was achieved by 
adding small amounts of Mg to the melt (<0.5% of the Bi 
molar amount). The carrier density of these samples was 
measured from the Hall and Shubnikov-de Haas effect to 
be 2.3x10^^ cm~^. The samples are cleaved in-situ un- 
der a base pressure < 1 x 10~^^ torr. We optically excite 
the sample with a linearly-polarized infrared laser pulse 
{hui = 1.5 eV, 50 fs) and subsequently probe the evo- 
lution of the transient populations in the occupied and 
unoccupied band structure by photoemitting electrons 
with an ultraviolet laser pulse {hu2 = 6.0 eV, 160 fs). 
The photoelectron kinetic energy and emission angle are 
resolved by a hemispherical electron analyzer. During 
measurement the sample temperature is maintained at 
70 K. Our laser setup is based on a Ti: Sapphire oscil- 
lator operating at 80 MHz repetition rate. In this par- 
ticular experiment, we compromise the ultraviolet pulse 
duration for a good spectral resolution of < 22 meV. The 
time resolution obtained by cross-correlation of hui and 
hiy2 is 163 fs. The absorbed infrared fiuence is computed 
to be 26 /iJ/cm^ using the measured reflectivity of 0.4 at 
45° incidence 0. 

Fig.[T](a) displays the transient photoemission intensity 
within the integration windows indicated in the subse- 
quent panels. A movie of the same data is included in the 
supplemental material [l3|. For reference, a schematic of 
the band structure is illustrated in Fig.[T](b). At negative 
delays — before the arrival of the infrared pump pulse 
— we probe the equilibrium electronic structure. Due to 
the p-type doping, the SS and BCB are unoccupied in 







I ■ I ■ I ■ I ■ I I ' I ' I ' I ' I I ' I ' I ' I ' I 

-10 -5 5 10 -10 -5 5 10 -10 -5 5 10 
Angle (deg) Angle (deg) Angle (deg) 



FIG. 1. (Color online) (a) Transient photoemission intensity 
within the integration windows indicated in the subsequent 
panels. Feature (1) corresponds to an IPS populated by /iz/2 
and probed by hui and thus appears before time-zero. Fea- 
ture (2) is a bulk state populated by a direct optical transition 
at time-zero. The BCB (3) and SS (4) are indirectly popu- 
lated by scattering from higher energy states, (b) Schematic 
of the electronic band structure for Bi2Se3. Grey, shaded re- 
gions represent the BVB and BCB. Lines represent the SS 
with spin-texture indicated, (c)-(g) trARPES spectra near 
the F-point for various pump-probe delays, (c) BVB before 
excitation, (d) Initial optical transition to high-lying bulk 
state and IPS. Note the separate energy axis for the IPS. (e) 
The SS and BCB are populated by scattering from higher- 
lying states, (f) Energy relaxation of the SS and BCB popu- 
lations are mostly complete, (g) A meta-stable population in 
the BCB, accompanied by a persistent population of the SS 
in the region energetically below it. A movie showi ng d ata at 
all delays is included in the supplemental material 13]. 



equilibrium, see Fig. [TJc). The first features to be popu- 
lated are displayed in Fig.[T]^d). Curiously, the parabolic 
dispersing feature (labeled 1) appears to be populated be- 
fore time-zero. Moreover, closer inspection reveals that it 
decays toward negative delays [IS] • We attribute this fea- 
ture to the first image potential state (IPS) [iJ, ll7H2QJ . 
which is an electron state bound in front of the metallic 
Bi2Se3 surface, and frequently observed in time-resolved 



photoemission experiments. Its decay toward negative 
delays indicates that it is populated by hu2 and probed by 
hui, hence its separate energy axis in the figure. We com- 
pute a binding energy Eips — E^ac = —0.77(3) eV, consis- 
tent with the IPS binding energy —0.85 eV expected for 
a perfectly metallic surface [13|, ll7|, |l8| . Importantly, the 
decay toward negative delay implies that the IPS relax- 
ation is decoupled from the SS and bulk band dynamics 
which are the focus of this Letter, and we can safely dis- 
regard it for the remaining discussion. 

The feature at E — Ep = 1.3 eV (labeled 2) decays with 
a time constant 70(20) fs and results from a direct optical 
transition pumped by hvi from the BVB edge to higher 
lying states in the bulk. Note that the lower- lying BCB 
and SS have negligible population around time-zero, and 
it takes ~ 700 fs for them to reach maximum population, 
see Fig.[T](a). This indicates that these bands are not di- 
rectly populated by hvi , but are populated indirectly by 
scattering from higher- lying states. The SS dispersion 
agrees with ARPES measurements in thermal equilib- 
rium, see Fig. [TJe) [y, |2l|, ensuring that we can discuss 
transient populations in a fixed electronic structure. The 
BCB does not have a sharp dispersion, consistent with 
its 3D bulk nature, but appears as a diffuse electron dis- 
tribution centered around the F-point. 

After 2 ps the SS and BCB populations have signifi- 
cantly decayed and energetically relaxed towards the bot- 
tom of their respective bands, see Fig. [Hf). The subse- 
quent dynamics is much slower and persists for > 10 ps. 
Obstructed from further decay due to the bandgap, the 
BCB electrons form a meta-stable population at the BCB 
edge. Intriguingly, this BCB population is accompanied 
by a persistent population in the SS, but only energeti- 
cally below the BCB edge, see Fig. [TJg). 

To understand the coupled dynamics of bulk and sur- 
face electrons in more detail, we proceed by analyzing the 
energetic distribution of electrons in the bulk. Fig. [2fa) 
presents energy distribution curves (EDCs), obtained by 
integrating the trARPES spectral intensity within a ±2° 
angular window around the F-point, for selected delays. 
We first focus on the energetic region of the EDCs asso- 
ciated with the BVB. By fitting each EDC with a Fermi- 
Dirac (FD) distribution, we extract an electronic temper- 
ature Tg as a function of delay, shown in Fig.[2fb) [13]. At 
negative delays there is a well-defined FD distribution at 
Ep in the BVB. After excitation, a FD with increased Tg 
is observed and attributed to scattering of photoexcited 
electrons with electrons in the cold Fermi sea [22|, |23|. 
The hot electron population subsequently cools by trans- 
ferring energy to the lattice [2J, |25|. Within 9 ps, Tg al- 
most returns to its equilibrium value, and the BVB elec- 
trons and lattice system can be regarded as equilibrated. 
We find that Tg decays exponentially with a time con- 
stant tbvb = 1.85(6) ps. Note that we confine our fit to 
delays after 2 ps to avoid non-analytic transients associ- 
ated with indirect filling processes. 
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FIG. 2. (Color online) (a) EDC at the F-point for various de- 
lays. The energetic regions associated with the BVB and BCB 
are specified, (b) The transient electronic temperatures ex- 
tracted by fitting Fermi-Dirac distributions to the EDCs. Af- 
ter 2 ps the BVB and BCB temperatures relax exponentially 
with time constants tbvb = 1.85(6) ps and tbcb = 1.67(5) ps, 
respectively. 



We now turn to the energetic region of the EDCs as- 
sociated with the BCB. The BCB is unpopulated before 
excitation, but within 300 fs a diffuse electron distribu- 
tion extends > 0.5 eV above Ep- The electrons are dis- 
tributed exponentially in energy, suggesting that intra- 
band scattering processes rapidly (within our 160 fs time 
resolution) redistribute the electrons within the band, 
taking them into a thermal distribution. The distribu- 
tion subsequently cools, again by transferring energy to 
the lattice. By fitting to a FD distribution we charac- 
terize the BCB population with an effective tempera- 
ture, shown in Fig. [2jb) [13]. Fitting Tg after 2 ps, as 
above, we find an exponential decay with time constant 
Tbcb = 1.67(5) ps. The near agreement of tbvb and 
Tbcb motivates us to identify their mean Te-p = 1.74(4) ps 
as the timescale for electrons to transfer energy to the lat- 
tice, and is thus a characteristic time for electron-phonon 
scattering in the bulk of Bi2Se3. We note that similar 
timescales for intra-band cooling have been observed in 
optical pump-probe measurements [llUl2U26|. 

After this energy-relaxation process, the system is 
characterized by a meta-stable BCB population. This 
long-lived BCB population is not unusual, and is in 
fact responsible for the phenomenon of photoconductiv- 
ity which is typically observed in semiconductors [27|. 
The longevity of this population is due to the fact that 
the electron-phonon scattering processes responsible for 
energy relaxation within the band are incapable of recom- 
bining electrons with the BVB. Here the highest phonon 
energy is 23 meV [28] , while the bandgap of 200 meV is an 
order of magnitude larger [21]. This slow recombination 
of bulk carriers in Bi2Se3 has previously been observed 
in optical measurements [26]. The simultaneous persis- 
tence of the metallic SS population, however, is surpris- 
ing. Photoconductivity of a metallic state is unexpected 
since the absence of a bandgap means there is no barrier 
to rapid recombination. We attribute this SS population 
to a continuous filling from the meta-stable BCB popu- 
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FIG. 3. (Color online) (a) Transient photoemission intensities 
within the integration windows indicated in the subsequent 
panel (b). Curve #10 corresponds to the integration window 
over the BCB edge, and is normalized to match the intensity 
of the SS intensities. The SS population decay is exponen- 
tial at energies above the BCB edge, but has a second, slower 
component below the BCB edge. This slow component decays 
with the same timescale of 5.95(2) ps as the BCB population, 
demonstrating that the SS is continuously filled by the slowly 
decaying BCB. (c-f) Schematic of the transitions and scatter- 
ing processes discussed in this Letter, including the (c) direct 
optical transition, (d) scattering into SS and BCB, (e) intra- 
band scattering of the BCB and SS, and (f) the BCB-to-SS 
scattering responsible for the persistent SS population. 



lation. The first evidence for this is Fig. [IJg), in which 
it is evident that the SS persists only energetically be- 
low the BCB edge. We find yet stronger evidence in the 
population dynamics. In Fig. [3l^a) we plot the transient 
photoemission intensity obtained by integrating within 
the windows indicated in Fig. [3l^b). Above the BCB, the 
SS decays with a single exponential. However, below the 
BCB a second slower component is observed, reflecting 
the presence of an additional filling channel. Moreover, 
the decay rate of this slower component matches that of 
the population decay at the BCB edge (see curves #6 
and #10). The agreement of these decay rates is direct 
evidence that the persistent SS population is due to con- 
tinuous filling from the meta-stable BCB population. We 
note that similar bulk-to-surface scattering behavior has 



been observed in Si, 
crossing Ep 



which however lacks a metahic SS 



2a, a 



We can deduce some properties of the transiently ex- 
cited populations relevant to transport applications. For 
bulk crystals in equilibrium, current is typically dom- 
inated by non-spin-polarized BCB carriers rather than 
the spin-polarized SS [9|,ll0|. In our experiment, the ratio 
of photoexcited SS carriers to BCB carriers is enhanced 
with respect to the equilibrium case due to the finite op- 
tical penetration depth of the excitation {a~^ = 50 nm) 
|l3| . Since the transient populations are generated only 
in this near-surface region, their transport properties are 
more comparable to those of thin films and nanoribbons, 
in which SS contributions are indeed observed [3l|, l32|. 
Thus, the photoexcited SS population could represent a 
channel through which to drive transient spin currents 
without overwhelming contributions from the bulk. 

An exponential fit to the BCB population decay (curve 
#10 in Fig.[3fa)) gives a population lifetime of 5.95(2) ps. 
Since the SS population is sustained by the BCB, it is rel- 
evant to investigate the mechanism for the BCB decay. 
We have argued that electron-phonon coupling cannot 
be responsible for recombination due to the low phonon 
energies. Direct recombination via photon emission is 
also a possibility, though the timescale for such a process 
is typically ;::^ 1 ns, and thus seems unlikely to be rele- 
vant on our measurement timescale [33|. The BCB decay 
could also be attributed to spatial diffusion of BCB elec- 
trons away from the probed volume of the sample. The 
diffusion coefficient D for Bi2Se3 has been measured to 
be 1.2 cm^/s [12]. Combining this with the absorption 
coefficient of the pump a = 0.02 nm~^ [135 , ^^ ^^^i" 
mate a characteristic diffusion time {a'^D)~^ = 21 ps 
[34] . This is rapid enough for diffusion to be a significant 
contribution to our measured lifetime. Characterization 
and control of spatial diffusion, for example by utilizing 
thin films with thickness less than the pump penetration 
depth, could therefore drastically extend the observed 
lifetime. 

In conclusion, we have performed trARPES measure- 
ments on p-type Bi2Se3 and elucidated the rich electron 
dynamics generated by optical excitation. We have iden- 
tified the first image potential state in front of the metal- 
lic surface, which may be of particular interest in this ma- 
terial due to the proposal that TIs host image magnetic 
monopoles [35] . We have also identified the direct transi- 
tions and subsequent relaxation processes induced by op- 
tical excitation, which we have schematically illustrated 
in Fig. [3l^c)-(f). Finally, we have shown that these tran- 
sitions culminate in a persistent non-equilibrium popula- 
tion of the spin-textured Dirac SS, which is attributed to 
continuous filling from a meta-stable population in the 
BCB. We note that excitation with any above-bandgap 
photon energy should lead to the same SS filling behav- 
ior, since inter- and intra-band scattering processes in- 
evitably lead to relaxation of carriers toward the BCB 



edge. This phenomenon could find a role in applications 
requiring ultrafast optical control of a spin-polarized sur- 
face conduction channel. 
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